A binary mixture of spinor atomic Bose-Einstein condensates 
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We study the ground state and classify its phase diagram for a mixture of two spin-1 condensates 
in the absence of external magnetic (B-) field according to atomic parameters for intra- and inter- 
species spin exchange coupling and singlet pairing interaction. Ignoring the inter-species singlet 
pairing interaction, the ground state phases are found analytically. Numerical approach of simulated 
annealing is adopted when the singlet pairing interaction is present. Our results on the phase diagram 
and the boundaries between phases allow for easy identifications of quantum phase transitions, that 
can be induced through the tuning of optical traps and atom numbers. They provide the first insight 
and guidance for several ongoing experiments on mixtures of spinor condensates. 
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The liberation of atomic hypcrfine spin from optical 
trapping opens the study of spin-dependent phenomena 
in atomic quantum gases [3, [2, 13] • In the simplest case 
of a spin-1 condensate, two quantum phases: polar and 
ferromagnetic, exist depending on the sign of atomic spin 
exchange interaction 0, 0, Q . The spin-2 system accom- 
panied by richer physics is significantly more complicated 
as in other higher spin systems 0, @, & [H3| ■ Recently 
several groups have initiated the study of spin-3 systems, 
of which the atomic 52 Cr is a viable experimental candi- 
date [111, [HJ ■ Parallel to the efforts on spinor condensates 
are active studies of condensate mixtures with more than 
one atomic species or state [3 Q H (3 (3 B El- 
Earlier works on double condensates rely on atoms with 
two almost identically trapped internal states [3, EH, 
while more recent experiments demonstrated mixed con- 
densates with tunable inter-species interactions [3 HH ■ 
Although the field is blossom with extensive studies on 
spinor condensates and mixtures of scalar condensates, 
few have touched the subject of mixtures of spinor con- 
densates [22|. 

This work concerns a binary mixture of spin-1 conden- 
sates, with the mixture of 23 Na and 87 Rb atoms being a 
special case [22[ . We study and identify the ground state 
phase diagram under the mean-field approximation and 
the single spatial mode approximation (SMA) for each of 
the two spinor condensates. Our theory can be extended 
to more general mixtures of higher spin condensates or 
mixtures with more than two constituents. 

The interaction between two different spin-1 atoms is 
described by the contact pseudo-potential V\ 2 {ri — f^) = 
(g^-Po + + Q^-PiWi ~ r 2 ), where g<£ 2 ] 2 = 

Airf^a^ 2 \/ fi, Oof 2 are s-wave scattering lengthes in the 
channel of total spin F to t = 0, 1, 2 respectively, /i is the 
reduced mass. Vo,i,2 is the corresponding projection op- 
erator. Using Fx ■ F 2 = Vi -V\ - 2V 0,11, we find 

V 12 (n - r 2 ) = (q + /3P! • F 2 + jVoJSin - r 2 ), (1) 
where a = (g[ 12) + g { 2 12) )/2, (3 = {-g[ 12) + g { 2 2) )/2, and 



7 = (2<7q — 3<7j +g 2 )/2. Between same species spin- 
1 atoms, the interaction takes the familiar form Vi j2 (ri — 
r 2 ) = (ai,2+/3i,2Fi ! 2-Fi j2 )(5(ri-r2), which respects the 
identical particle symmetry, thus does not include the 
odd symmetry term of total spin 1 projection 0, 

When atom numbers for both condensates are large, 
mean-field approximations can be satisfactorily applied. 
We further assume atomic interaction parameters are 
such that SMA 0, [23| for each spinor condensate holds, 
i.e., the mean- field mode functions are almost identi- 
cal for the three components of each spinor conden- 
sate. We thus take <l.,(r) = ^/Nl^if)^ and = 
y/N 2 4>(f)Q , respectively, for the mean fields of the two 
condensates. The mode functions ip(r) and 4>{r), de- 
termined presumably by the stronger density dependent 
interactions, are normalized according to J dr\ip(r)\ 2 = 
J dr\(t>(r)\ 2 = 1. Ni and N 2 denote total atom num- 
bers for the two condensates, distributed into the respec- 
tive spinor component with N^ 1 ' 2 ' (X^-^ 1 ^ = -^1,2)- 
The effective order parameters for the spinor fields are 

Cf 5 = \J /Nxe- 10 * and cf } = \J /N 2 e~ %Vi with 
0i and tfi real phase values. Neglecting spin-independent 
terms, the Hamiltonian for our model of a mixture of two 
spin-1 condensates takes the form 

H = ^NlCxfaft + \N 2 C 2 p 2 i 2 

+^iV 1 iV 2 C7 12 /3f 1 ■ f 2 + ijVii\TaCi 2 7|a-I 2 . ( 2 ) 

The interaction coefficients arc C\ = J dr\ip(r)\ 4 , C 2 = 
fdr\<j>(r)\ 4 , and C 12 = J dr\rp(r)\^{r)\\ f u = 

E /•(!)* rpi A 1 ) f _ V 1 A 2 )* pi A 2 ) „„J „ _ 
m,n Sm r mnSn , J2i — Z^m,n^ m r mnSn , 1:11111 b - — 

Y^, m {~ l) TO Cm C-ot- Whether the two condensates are 
miscible, i.e., completely overlap or not, is unimportant. 
As long as the SMA holds for each spinor condensate, our 
model retains its validity even for immiscible condensates 
with diminishing cross spin coupling oc C\ 2 . 

Given the interaction coefficients, the above Hamilto- 
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nian (J2J) describes a system of two interacting spins. It 
involves a total of six complex variables when treated 
semi-classically. The U(l) gauge symmetries for the 
two condensates give rise to two constraints: the con- 
servation of total atom numbers in each condensate. 
After a detailed calculation, three independent relative 
phases appear which we chose as: 7]± = 6\ + 8—± — 2#o, 
V2 =fi + if-i - 2(p , and 773 = 0_i + (pi - <f-i - 6\. 



To find the global ground state regardless of magneti- 
zation, we minimize Eq. with respect to \ N± , 
and the relative phases 771.2,3 defined above. Alterna- 
tively, we can minimize with respect to semi-classical 

spin variables. Both approaches give the same results, 
fi 2) 

provided that Q ' is replaced by the normalized value 

C| 1,2 V N^ 1 ' 2 ^ /Ni : 2 to take into account the con- 
straints explicitly. The multi-dimensional minimization 
is numeric ally carried out with the method of simulated 
annealing [2J] where analytical approach fails. Without 
loss of generality, we report the special case of JVi = iV 2 = 
N in the following. 

We first discuss the simple situation of 7 = 0. The 
results will serve as reference states for the complete 
ground state phase diagram. As the ground state is de- 
generate, it is convenient to restrict into the subspace of 
fi+ = fi- = (t = 1, 2), which leaves two unknowns: j\ z 
and / 2z ■ The ground state is analytically derived includ- 
ing explicit formulae for phase boundaries as functions of 
spin coupling parameters. Figure [1] illustrates the depen- 
dence of order parameters on the inter-species interaction 
Ci2/9 at fixed values of intra-species spin exchange inter- 
action parameter C\(3\ and Cifii-, which categorize the 
ground state into three cases. 

First, when the two spin-1 condensates are both fer- 
romagnetic, the ground state is simple, as shown in Fig. 
QJa). Both spins are polarized: ff = ff = 1. When 
Ci2/3 > 0, they are anti-aligncd as fi ■ f2 = — 1. When 
C12P < 0, they are aligned with fi • {2 = 1. 

Secondly when both condensates arc polar with 
C\(3\ > and C2P2 > 0, we assume C\(3i ^ C2P2 
without loss of generality. The ground state is found 
to exhibit several phases. If Ci 2 |/3| < 2v / CiAT^2, 
i.e., for weak intra-species spin exchange interaction, the 
two spins become essentially independent, giving rise to 
ff = 0, ff = 0, and fi • f 2 = 0. With increasing 
strength of C12 , one of the spins become polarized 
and the other is partially polarized for 2^Ci(3iC202 < 
Ci2\f3\ < 2max(Ci/3i, C-zfa)- The exact value for the 
resulting partial polarization is determined by the rel- 
ative strengths of the two intra-species spin-exchange 
interactions. For Ci/3i < C2P2, we have = 1, 
ff = (Ci2/3/2C 2 /3 2 ) 2 , and f 1 • f 2 = -C 12 f3 /2C 2 /3 2 . At 
stronger inter-species spin exchange interaction, when 
C 12 > 2max(Ci/3i, C2P2), the spin-spin coupling be- 
tween the two condensates causes all atoms to be polar- 
ized as illustrated in Fig. [Hb). 
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FIG. 1: (Color online). The dependence of ground state order 
parameters on at fixed values of Ci/3i and C2P2, for 

7 = 0. Blue solid lines, red dashed lines, and black dot-dashed 
lines denote respectively the mean field order parameters ff , 
ff , and fi • f2. The three subplots refer to the three cases of 
fixed intra-species spin exchange interactions (Ci/3i, 6*2/32)=: 
(a) (-1,-2); (b) (1,2); and (c) (-1,2). 



The third case as shown in Fig. [He), with one 
ferromagnetic (Cx/3i < 0) and one polar condensate 
{C2P2 > 0), is found to be similar to case (b) of two 
polar condensates. When Ci 2 |/?| < 2 max(Ci/?i, C2P2), 
the ground state gives ff = 1, ff = {C12P /2C2P2) , and 
fi • f2 = — Ci 2 /?/2C 2l 3 2 , which is the same as in case 
(b). When C l2 \P\ > 2max(Ci/?i, C 2 /3 2 ), we find ff = 1, 
ff = 1, and |fi • f 2 | = 1. We note that Luo et. al. [12], 
recently performed quantum calculations and studied the 
energy band structure of essentially the same model for 
this case using estimated inter-species atomic parameters 
for 87 Rb and 23 Na atoms HHg. 

Next, we discuss the general case of 7 7^ 0. Before con- 
sidering the full Hamiltonian , we illustrate the phase 
diagram for the special case of C\f3\ = C2P2 = in Fig. 
[2] Only two terms remain in the Hamiltonian, which for- 
mally resembles the Hamiltonian of a spin-2 condensate 
|| . The phases PP, FF, CC, and AA are characterized by 
the two order parameters (fi-f 2 = 0, |s_| = 1), (fi-f 2 = 1, 
= 1), (h ■ f 2 = 0, |s_| = 0), and (fi ■ f 2 = — 1, 
\s—\ = 1), respectively. The other two order parameters 
take fixed values (ff = 0, ff = 0), (ff = 1, ff = 1), 
and (ff = 1, ff = 1) respectively in the PP, FF, and AA 
phases, while for the CC phase, their values remain un- 
determined. Compared to a spin-2 condensate, we find 
the PP, FF, and CC phases are similar to the polar, fer- 
romagnetic, and cyclic phases respectively There exists 
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no counterpart for the AA phase in a spin-2 condensate. 
The two hypothetical (spin-1) sub-spins fi and f 2 making 
up the total spin-2 are constrained to be identical. In the 
case of a mixture as we discuss here, however, the two 
spins are different. The phase boundaries are specified 
by two lines with (3 = and 3/3 = 7 as in a spin-2 con- 
densate. The transitions between different phases are all 
first order in this case. 
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FIG. 2: (Color online). The ground state phase diagram for 
our model when only inter-species spin exchange and singlet 
pair pairing interactions are present. Red solid lines denote 
discontinuous phase transition boundaries. The PP, FF, CC, 
and AA phases are denoted by the values for the two order 
parameters (fi ■ f2,|s_|)=: (0,1), (1,0), (0,0), and (-1,1) re- 
spectively. 



The most general case of our model Hamiltonian @ is 
when all terms are present. As before for 7 = 0, we will 
restrict to respective cases specified by the fixed values 
for C\(3\ and C202- Using full numerical simulations, we 
find the dependence of ground state phases on the four 
parameters: Ci/?i, Cifii, C12P, and C127. The compe- 
titions among the first three determine the ground state 
as categorized in Fig. [TJ The 7-term encourages pairing 
two different type of atoms into singlets when 7 < 0; it 
minimizes |s_| instead when 7 > 0. If C12I7I is small 
enough compared to the first three parameters, then in- 
tuitively we expect the ground state structure of Fig. [T] 
will largely remain intact. However, our numerical results 
uncover remarkable differences between 7 > and 7 < 
even when its value is small. In the extreme circumstance 
when C 12 1/3| ^> 1 and C12I7I S> I, the ground state is 
fully determined by the respective ratio of (3 to 7, result- 
ing in the same phase diagram as for C\(5\ = C2P2 = 0, 
which is illustrated in Fig. O The resulting phases are 
categorized as before and illustrated in Fig. [3] 




FIG. 3: (Color online). The ground state phase diagram of 
our model at fixed values of C\(3\, C-ifii- Blue solid lines 
denote continuous phase transition boundaries. Red solid 
lines denote discontinuous phase transition boundaries be- 
tween two phases with fully determined order parameters. 
The green solid line denotes the discontinuous phase tran- 
sition boundary between the fully determined phase CC3 and 
the MM phase. The black dash-dotted lines correspond to 
C127 = 0, which serve as guides for the eye. The three sub- 
plots denote fixed intra-specie spin exchange interaction pa- 
rameters as in Fig. Q] (Ci/3i, C 2 /3 2 )=: (a) (-1-2); (b) (1,2); 
and (c) (-1,2). 



First, for two ferromagnetic condensates, correspond- 
ing to Fig. [IJa) of C127 = 0, we find two phases: AA 
and FF for C 12 /3 > and C 12 /3 < respectively. The 
phase transition between them are discontinuous. In- 
creasing C127 along the negative axis direction, we find 
the same two phases AA and FF if C12I7I remains small 
in value, although the phase boundary is found to shift 
with changing C12I7I. Further increasing the strength 
of singlet pairing interaction, a new phase MM emerges 
between the A A and FF phases, where all four order pa- 
rameters evolve continuously to border the AA and FF 
phases. Beyond a critical value of C12I7I, the PP phase 
occupies an interval section of the parameter C12/3. Next, 
we increase C127 from along the positive axis direction. 
In the beginning, a positive 7-tcrm decreases |s_|. No 
matter how small C12I7I is, the critical point (at 7 = 0) 
between the phases AA and FF diffuses into an interval, 



4 



and all four order parameters change continuously from 
the FF to the AA phase. When C127 is sufficiently large, 
as in the case of the opposite direction when C127 < 0, 
a special type of CC phase, which we call CC1 (f 2 = 0, 
ff = 1) arises over an interval of the parameter C12/?. 
A continuous transition region is found to surround the 
CC1 phase and borders the FF and AA phases. These 
results arc summarized in Fig. |3fa) , with the boundaries 
of continuous (discontinuous) phase transitions denoted 
by blue (red) solid lines. 

The case of two polar condensates (0 < Ci(3\ < C^fi-i) 
are shown in Fig. E^b). The C127 = line is seen to 
be partitioned into five intervals. The middle interval 
around Cyifi = is the region where the two spins are 
completely independent as discussed before. Addition- 
ally, two continuous changing intervals belong to the MM 
phase. The remaining two intervals are the A A and FF 
phases, the same as in the first case of Fig. EJa). The 
inclusion of the fourth competing interaction 7-term, we 
find a new phase emerges on each side of the C127 = 
line irrespective of the value of the 7-term. This confirms 
the middle interval is really a boundary between the PP 
phase and a second special type of CC phase, denoted 
as CC2 (f 2 = 0, ff = 0). The transition across is found 
to be discontinuous. When C127 < 0, the PP phase is 
found to dominate. In the opposite case, the CC2 phase 
wins. 

The third case corresponds to a mixture of one ferro- 
magnetic (C1P1 < 0) and one polar (C2P2 > 0) conden- 
sate. When C127 = as discussed before, there exist two 
phases AA and FF in both sides of the C12P axis, plus a 
middle continuous transition region. When C127 7^ 0, we 
first decrease C127 along the negative axis from 0. Before 
reaching a critical value, the ground state phase is found 
to be essentially the same as for C127 = 0. Beyond this 
critical value, a PP phase emerges as in the first case of 
Fig. E^a). Increasing C127 7^ along the positive axis, 
the results is found to change completely. The critical 
value disappears. A third type of CC phase, labeled as 
CC3 (ff = 1, ff = 0), emerges in the middle region. The 
transition between the CC3 and the middle transition 
region MM phase is discontinuous on both boundaries 
along the axis of C12P, denoted as green solid lines in 
Fig. E|c). This discontinuous phase transition boundary 
is distinct from that of the red solid lines, which denote 
boundaries between two phases with fixed order parame- 
ters and the changes for all four order parameters across 
the boundaries are exactly known. 

Finally we hope to point out that although inter- 
species interaction parameters are fixed, the interaction 
coefficients in the model wc consider such as C12P and 
C127 remain tunable by controlling optical trapping po- 
tentials or adjusting atom numbers. In this study, we 
consider the special case of two equally populated con- 
densates inside fixed traps. In reality the effect of differ- 
ent atom numbers can be used to adjust the four interac- 
tion coefficients. We can also use different optical traps 
for the two species of atoms, which can further enhance 



the tunability of the different ranges and ratios of the 
values for C\, C%, and C\2- 
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FIG. 4: (Color online). The dependence of the ground state 
order parameters on the trap frequency ui^a,, at a fixed u-n.b = 
2ir x 50Hz. Blue squares, red circles, and black diamonds 
denote respectively for the three order parameter fj 2 , f|, and 
fi • f2 . Lines are guides for the eye obtained from connecting 
neighbouring data points. Dot-dashed line specifies the phase 
boundary between the MM and the FF phase. 



As a concrete example, we consider here a realizable 
experiment to demonstrate a quantum phase transition 
from the MM phase to the FF phase in a binary mix- 
ture of 87 Rb (species one) and 23 Na (species two) spin-I 
condensates. The tuning is achieved through the con- 
trol of the trapping frequency for Na atoms. We assume 
5 x 10 5 87 Rb atoms and 10 5 23 Na atoms are confined 
harmonically V a {r) = ^M a uj a r 2 with = 2tt x 50Hz, 
and WNa tunable from 2ir x 30 to 2tt x 80Hz. The inter- 
species atomic parameters between 87 Rb and 23 Na are 
unknown. We therefore resort to the simple approach 
of degenerate internal states [25|, HH, [27} as the low en- 
ergy atomic interactions can be mostly attributed to the 
ground state configurations of the two valence electrons. 
The inter-species scattering lengths for singlet and triplet 
electronic states are approximately determined already, 
given by as = 109ao and Or = 70ao [25l. l2r|, where ao is 
the Bohr radius. The inter-species interactions between 
87 Rb and 23 Na atoms are then parameterized by three 
scattering lengths, each being a linear combination of as 
and ar weighted by the appropriate 9j coefficients for the 
total combined spins of F to t = 0, 1, and 2 I22| . Within 
this approximation, it's found coincidently [22j that 7, 
the parameter for inter-species singlet pairing interaction 
is equal to zero. The ground states are found through 
imaginary-time propagation of the corresponding cou- 
pled Gross-Pitaevskii equations. To compare with our 
former results within SMA, we redefine the three order 
parameters as f 2 = / dr(ip*Fijipj) 2 / J dr(ip*ipi) 2 , ff = 
/ dfitiFijfa) 2 / J df{<f>*<Pif, and fi-f 2 = / dr^F^)- 
(ftFij^j)/ 1 dr(if)*ipi) ((/)*&), where ipi and fa (i = 
1,0,-1) are component wave functions respectively for 
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87 Rb and 23 Na condensates. When the SMA is satisfied, 
they reduce to that as defined in Eq. @, and the values 
of the order parameters become the same as in Fig. QJc) . 
This is indeed the case for most of our numerical calcula- 
tions, especially in the FF phase, with all atoms fully spin 
polarized, or when two condensates are uncorrelated. 

In Fig. 31 the three order parameters for the ground 
state are demonstrated under the changing trapping fre- 
quency of 23 Na atoms. With increasing value of WN aj the 
ferromagnetic spin-spin coupling between 87 Rb and 23 Na 
atoms will induce the spin polarization of 23 Na atoms, 
from initially unpolarized to a full polarization. Com- 
pared to the phase diagram of a mixture with one ferro- 
magnetic and one polar spin-1 condensate as in Fig. [3]Jc), 
it's clear that by tuning the trap frequency for the 23 Na 
condensate, we induce a quantum phase transition from 
the MM phase to the FF phase in the mixture. 

To conclude, we study the ground state of a binary 
mixture of two spin-1 condensates in the absence of B- 
field, using both analytical approaches and the numeri- 



cal method of simulated annealing. Neglecting the sin- 
glet pairing interaction between the two species of atoms 
(7 = 0), we present analytical formulae for all order pa- 
rameters of the ground state phase diagram including all 
phase boundaries and the development of order param- 
eters. When 7 ^ 0, we find the phase diagram for the 
ground state at three groups of specific values for the 
intra-species spin coupling interactions. We suggest it is 
possible to tunc to different phases, or induce quantum 
phase transitions by changing optical traps and/or atom 
numbers in the two condensates respectively. Finally, we 
confirm that the mean field theory we adopt remains an 
excellent approximation as all our results have been re- 
peated numerically with the full quantum approach of 
exact diagonalization at small atom numbers. 

This work is supported by US NSF, NSF of China un- 
der Grant 10640420151 and 10774095, and NKBRSF of 
China under Grants 2006CB921206, 2006AA06Z104 and 
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